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In this paper, we present the theory for cal- 
culating Raman line shapes as functions of 
the Fermi energy and finite temperatures 
in zinc blende, n-type GaAs for donor den- 
sities between 10 16 cm -3 and 10 19 cm -3 . 
Compared to other theories, this theory is 
unique in two respects: 1) the many-body 
effects are treated self-consistently and 2) 
the theory is valid at room temperature for 
arbitrary values of the ratio R = (Q 2 /a), 
where Q is the magnitude of the normal- 
ized wave vector and a is the normalized 
frequency used in the Raman measure- 
ments. These calculations solve the charge 
neutrality equation self-consistently for a 
two-band model of GaAs at 300 K that 
includes the effects of high carrier concen- 
trations and dopant densities on the per- 
turbed densities of states used to calculate 
the Fermi energy as a function of tempera- 
ture. The results are then applied to obtain 
the carrier concentrations from Fermi 
energies in the context of line shapes in 
Raman spectra due to the coupling 



between longitudinal optical phonons and 
plasmons. Raman measurements have 
been proposed as a non-destructive 
method for wafer acceptance tests of carri- 
er density in semiconductor epilayers. The 
interpretation of Raman spectra to deter- 
mine the majority electron density in n- 
type semiconductors requires an interdisci- 
plinary effort involving experiments, theo- 
ry, and computer-based simulations and 
visualizations of the theoretical calcula- 
tions. 
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1. Introduction 

The carrier concentration is a key figure of merit 
associated with a go/no-go decision for determining 
whether a wafer or an epitaxial layer meets specifica- 
tions and should undergo further processing. 
Technology roadmaps from the microelectronics and 
nanomaterials industries [1-3] call for non-destructive 
and fast-turn around methods to measure transport 
properties such as carrier concentrations in semicon- 
ductor wafers and epitaxial layers. Non-destructive 
measurements are economically more significant for 
III-V compound semiconductor wafers with epitaxial 



layers than for Si-based wafers with epitaxial layers 
because the regions of wafers used for making contacts 
cannot be used for product. Contacting measurement 
methods may be acceptable for Si wafers, but such 
destructive methods are much less acceptable for III-V 
compound semiconductor wafers with epitaxial layers 
because Si wafers are much less expensive per unit area 
than compound semiconductor wafers. 

Raman spectroscopy has been proposed as one pos- 
sible way to measure carrier concentrations non- 
destructively The shape of the Raman spectral lines 
due to the longitudinal-optical phonons interacting with 
the plasmon collective modes of the electron gas (so 
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called coupled longitudinal optical (LO) phonon-plas- 
mon modes) provides information on transport proper- 
ties of the electron gas in polar semiconductors [4-8]. 
The frequencies CO of the coupled modes are propor- 
tional to carrier concentrations, and the peak widths Aa> 
of the coupled modes are proportional to the scattering 
rates due to electron-phonon interactions. The qualita- 
tive determination of carrier concentration and mobili- 
ty, which is inversely proportional to the scattering rate, 
from Raman spectra is reasonably straightforward 
based on these proportionalities. The quantitative deter- 
mination of carrier concentrations and mobilities 
requires more sophisticated modeling of the spectra. 
Many of these higher level models involve fitting the 
spectra with the Fermi energy as a parameter and then 
determining the carrier concentration from knowing the 
fitted Fermi energy [6]. 

Most interpretations of Raman measurements on 
compound semiconductors such as GaAs require phys- 
ical models and associated input parameters that 
describe how carrier densities vary with dopant concen- 
trations and Fermi energies. In this paper, we develop 
the theory for extracting electron densities from Raman 
measurements of n-type GaAs at room temperature. We 
introduce two main classes of models that relate carrier 
concentrations to the Fermi energy for a given temper- 
ature and donor dopant density: 

a) Bandgap narrowing (BGN) models based on two 
equivalent bands at the T point in the first Brilluoin 
zone, and 

b) Parabolic densities of states (PDOS) models with 
and without a quartic term in the electron energy dis- 
persion E cT (k) for the T conduction sub-band, where k 
is the wave vector. 

The BGN models include many-body quantum 
effects and bandgap narrowing. The many-body quan- 
tum effects describe the electron-electron, electron- 
hole, and dopant ion-carrier interactions. The results 
are unique in two respects: 1) the many-body effects 
are treated self-consistently and 2) the theory is valid at 
room temperature for arbitrary values of the ratio R = 
(Q 2 /cc), where Q is the magnitude of the normalized 
wave vector and a is the normalized frequency used in 
the Raman measurements. Other reported work either 
is valid at low temperatures near K for arbitrary ratios 
(Q 2 /cc) [4] or is valid at arbitrary temperatures for ratios 
R typically R « 1 [5-7]. 



2. Theory 

Because the Fermi energy is one of the variables for 
calculating line shapes in the Raman spectra from semi- 
conductors, we first present in Sec. 2.1 the theory for 
the BGN and PDOS models by which numerical meth- 
ods give closed-form analytic expressions that relate 
carrier concentrations to Fermi energy E F , temperature 
r, and donor dopant densities N D . Then, within the con- 
text of these BGN and PDOS models, we develop in 
Sec. 2.2 the theory for calculating the line shapes 
observed in Raman spectra. 

2.1 Fermi Energy and Electron Density 

The temperature and donor dopant density are the 
independent input parameters for evaluating the Raman 
line shape factor given in Eq. (19). The Fermi energy 
depends in turn on the dopant density, in this case the 
donor density 7V D and T. For given 7V D and T, we pres- 
ent here models to compute the Fermi energy with full 
Fermi-Dirac statistics for the carriers at finite tempera- 
ture. By using these models, we calculate the majority 
electron density as a function of the Fermi energy in 
zinc blende, n-type GaAs for donor densities between 
10 16 cm -3 and 10 19 cm -3 . These calculations solve the 
charge neutrality equation self-consistently in terms of 
two main classes of models: the bandgap narrowing 
(BGN) model and three different parabolic densities of 
states (PDOS) models. 

Some researchers propose that Raman spectra may 
be a way to determine temperatures with spatial resolu- 
tions across wafers on the order of micrometers when 
the Fermi energy and dopant density are known from 
other independent measurements. Comparing the 
results given by the BGN and PDOS models, which are 
described below, also provides predictions on the sen- 
sitivity of the Fermi energy and Raman spectrum to 
variations in temperature. 

2.1.1 BGN Model 

The bandgap narrowing (BGN) model is a two-band 
model with one equivalent conduction band and one 
equivalent valence band at the T point in Brilluoin 
space. The BGN model is related to earlier work on n- 
type GaAs [9] and includes modifications to the densi- 
ties of states due to high concentrations of dopants, 
bandgap narrowing, and many-body effects associated 
with carrier-carrier interactions (carrier-carrier 
exchange and correlation). This BGN model is fully 
self-consistent and uses the Klauder self-energy (fifth 
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level of approximation) to calculate the distorted-per- 
turbed densities of states for the carriers. This BGN 
model has the following main features: 

1 ) Many-body quantum effects of carrier-carrier inter- 
actions and carrier-dopant ion interactions, bandgap 
narrowing, and distorted-perturbed densities of states 
for the carriers [9], 

2) Iterative and self-consistent solutions of the coupled 
charge neutrality equation and Klauder's fifth level of 
approximation for the renormalized self-energy propa- 
gator from which the distorted-perturbed densities of 
states are calculated, 

3) Full Fermi-Dirac statistics for the carriers at finite 
temperature, and 

4) Statistical analyses to give closed-form analytic 
expressions from very large, calculated data sets for 
carrier densities as functions of the Fermi energy 
Tables 1 and 2 contain the input parameters for the 
BGN calculations in Ref. [9] and for the BGN model 
given here. 

The electron n and hole h concentrations in units of 
cm -3 at thermal equilibrium are given, respectively, by 

n = J" /„ (E) Pc (E)dE and h = \ "[1 - f ( E)] p v ( E)d E, 
J- J- (i) 



where p c (E) and p v (E) are, respectively, the electron 
density of states for the equivalent conduction band and 
the hole density of states for the equivalent valence 
band [9], where f (E) = {1 + exp[(£ - E v )lk B T\Y l is the 
Fermi-Dirac distribution function. The calculations 
incorporate the Thomas-Fermi expression for the 
screening radius, 



\ne 



ee n 



> df (E ) 
dE 



[p c (E)-p v (E)]d% (2) 



and the charge neutrality condition 

N 1 =n-h, 



(3) 



to compute self-consistently the Fermi energy E F and 
the screening radius r s for given values of the ionized 
dopant concentration N Y and temperature T. The static 
dielectric constant is e and the permittivity of free space 
is £q. The ionized dopant concentration is positive for n- 
type material (donor ions) and negative for p-type 
material (acceptor ions). The results reported here are 
for n-type material. For the case discussed here, Nj = 
N D . The results for the screening radius r s are not 
reported here because they are not needed to extract 
carrier concentrations from Raman scattering measure- 
ments. 



Table 1. Fundamental Constants 



Parameter 


Symbol 


Value 


Units 


Planck's constant 


n 


6.5836 x 10" 16 


eV-s 


Boltzmann's constant 


k B 


8.6174 x 10" 5 


eV/K 


electron rest mass 


m 


9.1072 x 10" 28 


g 


electronic charge 


e 


^.802 x 10" 10 


esu 


Bohr radius 


a B 
e 2 /2a B 


0.5291 x 10" 8 


cm 


energy associated with 1 Rydberg 


13.6 


eV 


speed of light 


c 


2.9979 x 10 10 


cm/s 


wave length associated with 1 eV 


M 


1.2396x10^ 


cm 


wave number associated with 1 eV 


W 


8.0668 x 10 3 


cm -1 


dielectric constant in vacuum 


£ 


8.854 x l(T 12 


F/m 



Table 2. BGN model input parameters for intrinsic zinc blende GaAs at 300 K. The energies of the extrema of the con- 
duction and valence sub-bands are referenced to the bottom of the conduction sub-band at the r symmetry point in the 
Brillouin zone of the reciprocal lattice space. The mass of the free electron is m Q . These GaAs data are from Ref. [12]. 



Parameter 



Symbol 



Value 



Units 



bandgap 

effective mass for conduction band (2-band model) density of states 

effective mass for valence band (2-band model) density of states 

number of equivalent conduction bands 

number of equivalent valence bands 



£ G = |-£vrl 

m c 

m y 

N c 



1.424 

0.067 

0.572 

1 

1 



eV 
m 
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2.1.2 PDOS Models 

The PDOS models use parabolic densities of states 
for all equivalent bands and sub-bands. Unlike the 
BGN model in Ref. [9], PDOS models such as the four- 
band PDOS model for GaSb in Ref. [10] do not include 
modifications to the densities of states due to many- 
body effects and high concentrations of dopants and 
carriers because of computational limitations associat- 
ed with treating a four-band model in the context of the 
Klauder self-energy method (fifth level of approxima- 
tion). 

The zero of energy is at the minimum energy value 
(bottom) of the conduction T sub-band, E cT0 = 0.0. The 
bottoms of the conduction L and X sub-bands are, 
respectively at E cL and E cX . The maximum energy value 
(top) of the degenerate valence T sub-band is -E G , 
where E G is the intrinsic bandgap of GaAs. The split-off 
valence sub-band at r due to spin-orbit coupling is neg- 
lected. The probabilities for typical holes in equilibri- 
um to occupy appreciably these states in the split-off 
valence sub-band at r are very low. This means that the 
Fermi energies should be sufficiently above the valence 
sub-band maximum at T, Placing exact limits on the 
Fermi energies for which the PDOS models are valid 
would be tenuous, because knowledge of how the vari- 
ous sub-bands move relative to one another due to the 
dopant concentrations considered here and due to many 



body effects is not adequate. Table 3 lists the input 
parameters for the PDOS models. 

The heavy hole mass m^ and light hole mass m^ for 
the two degenerate sub-bands at the top of the valence 
band are combined to give an effective mass 



/ 3/2 . 3/2\2/3 

W vK=Ohh +"hk 



2 f 



(4) 



for the valence topmost sub-band, which becomes the 
equivalent valence band with a hole energy dispersion 
given by E vT (k) **-E G - (ft 2 I^/2m vT m ). 

The general expression for the parabolic densities of 
states for electrons and holes per band extrema and per 
spin direction is given by 



P {E): 



NAnVjE 



($K 3 )(h 2 /2m*m y 



(5) 



where N e is the number of equivalent ellipsoids in the 
first Brillouin zone, the volume of the unit cell is V= 
a L 3 , a L is the lattice constant, m* is one of the effective 
masses listed in Tables 2 and 3 for the appropriate band 
extrema, and m is the free electron mass. 

We sub-divide the PDOS models into the PDOS2, 
PDOS2NPG, and PDOS4 models. All three PDOS 
models include the equivalent valence band described 
by Eq. (4). 



Table 3. PDOS model input parameters for intrinsic zinc blende GaAs at 300 K. The energies of the extrema of the conduction and valence sub- 
bands are referenced to the bottom of the conduction sub-band at the r symmetry point in the Brillouin zone of the reciprocal lattice space. The 
mass of the free electron is m . These GaAs data are from Refs. [11,12]. 



Parameter 



Symbol 



Value 



Units 



bottom of the conduction L sub-band 

bottom of the conduction X sub-band 

top of the degenerate valence r sub-band 

spin-orbit splitting 

top of the split-off (spin-orbit splitting) valence r sub-band 

effective mass of conduction r sub-band 

non-parabolicity factor (quartic term prefactor) for conduction r sub-band 

transverse L sub-band mass 

longitudinal L sub-band mass 

effective mass of conduction L sub-band 

transverse X sub-band mass 

longitudinal X sub-band mass 

effective mass of conduction X sub-band 

light hole mass of valence r sub-band 

heavy hole mass of valence r sub-band 

effective mass of valence r sub-band 

splitoff band mass of the valence sub-band at r 

number of equivalent conduction L sub-bands 

number of equivalent conduction X sub-bands 



E G = \-E yr \ 


1.424 


eV 


EcL 


0.29 


eV 


E vX 


0.48 


eV 


~E v r 


1.424 


eV 


-E*, 


0.34 


eV 


~ E soT = ~ E vF - E$o 


1.764 


eV 


m cV 


0.063 


m 


5 


0.824 




«tL 


0.075 


m 


, 2,1/3 
m cL=OlLf%J 


1.9 
0.222 


m 
m 


™tX 


0.19 


m 


m cX=OlX m tx) 


1.9 
0.409 


m 
m 


% 


0.082 


m 


m W 


0.51 


m 


"Vr 


0.53 


m 


™so 


0.15 


m 


*cl 


4 




A^cx 


3 
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PDOS2 Model 

The PDOS2 model uses one equivalent conduction 
band and one equivalent valence band at the T symme- 
try point in the Brilluoin space for the integrals that 
appear in Eqs. (1) to (3). The electron energy dispersion 
for the equivalent conduction band is E cT (k)~E cT0 + 
(n 2 ^/2m c m ). 

PDOS2NPG Model 

The PDOS2NPG model is a two-band model with 
one equivalent conduction band and one equivalent 
valence band at the T point in Brilluoin space. It has no 
bandgap narrowing, but it includes the non-parabolici- 
ty for the electron energy dispersion in the equivalent 
conduction band at T. According to Ref. [11], we may 
include non-quadratic \k\ l terms in the electron energy 
dispersion E cT (k) for the conduction T sub-band in 
GaAs when k is small, namely, 

(6) 

where £ is the non-parabolicity factor. We use the Kane 
three level kp model [11], which does not include the 
conduction sub-bands at L and X, to include quartic 
terms in E(k) with 1 = 4. 

PDOS4 Model 

The PDOS4 model has three conduction sub-bands 
at the respective r, L, and X symmetry points in the 
Brilluoin space and one equivalent valence band at the 
T symmetry point. For the PDOS4, we modify here the 
PDOS model for GaSb in Ref. [10] so that it is valid for 
GaAs. It uses the parabolic electron energy dispersion 
E cT (k) for the conduction T sub-band in GaAs when k 
is small, namely, E cV (k) ~ E cro + {f^]^l2m Q m^. The 
general expression for the temperature dependence of 
conduction sub-band minima relative to the top of the 
valence band at T is given by [12], 



E i =E m -[A i T 2 /(T + B i )] 



(7) 



in units of eV, where / = T, L, or X. The values for the 
coefficients E lV ,A i9 and B l are listed in Table 4. Because 
8 permutations of the wave vector in the (1 11) direction 
exist, there are 8 L sub-band ellipsoids with centers 
located near the boundary of the first Brillouin zone. 
Also, because 6 permutations of the wave vector in the 
(100) direction exist, there are 6 X sub-band ellipsoids 



Table 4. Coefficients for the temperature dependence of the conduc- 
tion band extrema that are used in Eq. (15). These data are from Ref. 
[12]. 



Parameter 



Symbol Value 



Units 



r sub-band coefficients 



L sub-band coefficients 



X sub-band coefficients 



E T0 


1.519 
5.405 x 1(T* 


eV 
eV/K 


B r 


204 


K 


£ L0 


1.815 
6.05 x 10^ 


eV 
eV/K 


B L 


204 


K 


E xo 


1.981 


eV 


^ 


4.60 x 10^ 


eV/K 


By 


204 


K 



with centers located near the boundary of the first 
Brillouin zone. Since about half of each ellipsoid is in 
the neighboring zone, the number of equivalent sub- 
bands N cL for the E cL is four, and the number of equiv- 
alent sub-bands N cX for the X sub-band is three. 

In terms of the four-band PDOS4 model for room 
temperature n-type GaAs, the total density of states 
p c {E) for the majority carrier electrons in n-type GaAs 
then becomes 



p c (E) = p cr (E) + p cL (E) + PcK (E), 



(8) 



where p cT (E), p cL (E), and p cX (£) are the sub-band den- 
sities of states for the conduction T, L, and X sub-bands 
with effective masses of m cr , m cL , and m cX , respective- 
ly. The density of states for the minority carrier holes is 

p w {E) = p vY {E) (9) 

with an effective mass ofm vr . 

2.2 Dielectric Response Function 

The longitudinal optical (LO) phonons and plasmons 
interact in polar semiconductors such as GaAs to form 
LO phonon-plasmon modes. The theoretical line shape 
function L A (q,G)) of the Raman spectrum due to longi- 
tudinal optical (LO) phonon-plasmon coupled modes is 
then given for the configuration to which this theory 
will be applied by [8] 



L A (9,0)) = ^-e 



-h(o/k B T 



vV ft^-ay 2 II i_l 

> ^ {] -cr J ,n | E{q,6t)f 



fioi 
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where q is the scattering wave vector, (O is the Raman 
angular frequency shift, fi is the Planck constant, k B is 
the Boltzmann constant, T is the temperature in Kelvin, 
O) = G>ro(l + Qh) 1/2 is a parameter with the dimensions 
of angular frequency C FH is the dimensionless Faust- 
Henry coefficient which includes the LO/transverse 
optical (TO) phonon scattering ratio, and £%, is the TO 
phonon angular frequency Table 5 contains values for 
the parameters in Eq. (10). The total longitudinal 
dielectric response function £(#,&)) in the random 
phase approximation (RPA) is described by 

£(q,co) = \ + 4KX VE +4KX L (co)+4KxM c °l GO 

where the dielectric susceptibility Xve * s the contribu- 
tion from valence electrons, X h (G)) is the contribution 
from the polar lattice phonons, and Xefa *) is the con- 
tribution from the conduction electrons. The high fre- 
quency dielectric constant £ M is defined to be £„ = 1 + 
4/t^ye. Equation (11) then becomes for a binary semi- 
conductor, 



£(q,G>) = £» +47# L (fi>) +4?%(« <»)• 



(12) 



trie response function 47T£ e (</,6)) that describes light 
scattering by the conduction electrons in doped semi- 
conductors. We introduce the collision relaxation time 
T that describes the losses associated with electron- 
phonon and electron-dopant interactions and the corre- 
sponding angular collision frequency r= T" 1 . Mermin 
[14] showed that replacing co with (D+ir in the 
Lindhard expression for 47T£ e (</,6)) fails to conserve 
the number of local electrons and therefore is not the 
correct way to include collision broadening in 
47T£ e (</,6)). Instead, Mermin assumed that within the 
framework of a relaxation time approximation the elec- 
tron-phonon and electron-dopant interactions relax the 
electron density matrix to a local equilibrium density 
matrix [4, 14]. He then obtained the following 
Lindhard-Mermin relation [14] for ;£,(</, 6)), 

®Xs (q,0)+ir%*(q,a> + ir) 

where the Lindhard expression XeiQ^ + z ^") f° r e l ec ~ 
trons occupying states in a single equivalent conduction 
band is given by, 



The contribution of the polar lattice is given by [8] 



4;# L (©)=£„ 






(13) 



when phonon damping may be neglected and where 
0>l is the LO phonon angular frequency 

Within the context of the RPA, the Lindhard expres- 
sion [13] gives the electronic contribution to the dielec- 



4^ e °(«,0)+ir) = - rT J/ (£) 

71 q 



Ti 2 q 2 ft 2 qk _ / .„ 



2m c m m c m Q 



( fi 2 q 2 h 2 qk ., ._. 

+ — — +h((o+ir) 

I 2m c m Q m c m Q 



d 3 k, 



(15) 



Table 5. Dielectric response function input parameters for intrinsic zinc blende GaAs at 300 K. The energies 
of the extrema of the conduction and valence sub-bands are referenced to the bottom of the conduction sub- 
band at the r symmetry point in the Brillouin zone of the reciprocal lattice space. The mass of the free elec- 
tron is m . These GaAs data are from Ref. [14]. 



Parameter 



Symbol 



Value 



Units 



lattice constant 

static dielectric constant 

high frequency dielectric constant 

longitudinal optical (LO) phonon energy 

transverse optical (TO) phonon energy 

Energy associated with the angular collision frequency r 
due to electron-phonon and electron-dopant ion interactions 
Faust-Henry coefficient 

effective mass for the single equivalent conduction 
band density of states in Eq. (6) 



«L 


5.65 x 10" 8 


cm 


3) 


13.1 




£ M 


10.9 




hwto 


0.0353 


eV 




285 


cm 


hG>j 


0.0332 


eV 




268 


cm 


hr 


-0.0124 


eV 




-100 


cm 


Qh 


-0.4 




m c 


0.067 


m o 
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where fo(E) = {1 + exp[(£-£' F )/£ B 7]}~ 1 is the Fermi- 
Dirac distribution function, E F is the Fermi energy in 
eV, and T is the temperature in Kelvin. 

The integrand in Eq. (15) is independent of the 
azimuthal angle <p so that, 



The theoretical line shape function L A (q,cd) of the 
Raman spectrum then becomes 

L A (q 9 (0) 



1111UH.11C41 C411 6 1V Y ™ L11C4L, / 2 2 Y [ / N 1 



where jx = cos#, djti = -sin#, and qk = qkjx. To simplify 
Eq. (15) further, we introduce the following dimension- 
less normalized quantities: Q = qa B , K = ka B , K(E) = 

^(2m m c a B 2 /h 2 )E , and Q 2 = {h(w+ zT)m c /13.6 eV} 
where (e 2 /2a B ) = (fi 2 /2m a B 2 ) = 13.6 eV. Equation (15) 
then becomes with these definitions the following 
expression, 

J 1 dju(l+exp[{(13.6 eVK 2 lm c ) -E V }1 k B T\) 1 

x[(e 2 +2e^-Q 2 ) _1 +(e 2 -2e^+^ 2 )" 1 } ( 16 ) 

Performing analytically the angular integration over jti 
gives the complex results 

4?% (?,G>+ir) 

=^J o "^d^(l+exp[{(13.6 €MK 2 im c ) -E F }/ ^T])" 1 

7TC/ 



In 






+ ln 



c+r 



t_ 2 +7 2 



+ /(?? + -??_) + /(0- + -C7_) , 

(17) 



where p ± = g 2 ± 2QK- a, t ± = £> 2 ± 2g^ + a, 
ifc = -sign(p ± )arctan(7/p ± ) - (tt /2)[1 - sign(p ± )], 
a ± = +sign(T ± )arctan(y/r ± ) + (tt /2)[1 - sign(r ± )], 
Q 2 = a+i% a=(hom Q l\3.6 eV), y=(/?r#w c /13.6 eV), 
and -n < arctan(y/) < ^. 

When the temperature T=0 9 the Fermi function 
f (E) = (l +exp[{(13.6 eV K 2 im c )-E v }ik B T\y l is the 
unit step function. Then, the integral in Eq. (17) 
J ^ 2 dAT(...) is j F ^ 2 dAT(...) ; and an analytic eva luation 
is possible [4], where K F = yJm c E ¥ /13.6 eV . But, 
when T is room temperature, analytic evaluations are 
not possible. Researchers approximated the integral 
J K 2 dK(...) at finite temperatures by expanding its inte- 
grand in terms of the ratio R = (Q 2 /cc) for either very 
small or large R. However, such expansions do not cap- 
ture all of the spectral information in Eq. (17). 



(18) 



where Inn 



e(q,(Q) 
e T {q,m) = Re{e(q, ©)}= £, 
+4?% ; («,®+*n, and 



\e 2 {q,a>)+^{q, 

1 «ho 



ft>) 



-°ho 

-ft) 2 



e i (q,to)=Tm{e(q,to)}=4n% a (%to + ir)- 

The Lindhard-Mermin relation Eq. (14) then gives 
expressions for ^ re (?,ft)+ /I") = Re{x e {q,(0+ ir)} and 
^; e (?,ft>+ ir) = \m{x e (q,(Q+ ir)} as functions of 
Xj{q,a* + ir) = Re{ Xe \q,CO + ir)} and ^(g,(0 + ir) 
= Im{^ e °(^r,ft)+/F)}. Even though the factors in Eqs. 
(10) and (18) give the contributions to the Raman line 
shape for each scattering mechanism, an examination 
of only the factor 



Im< 



s(q,co)\ J 



£/(?,«>) 



e t 2 (q,co)+e i 2 (q,co) 



i Ko 2 -C0 2 ) 2 47tX K (q,(Q+in 
I D 2 (q,(0+ir) 



(19) 



gives the spectral information for extracting the Fermi 
energy. The denominator in Eq. (19) is given by the 
expression, 

D 2 (q,co+ir) = eJ((o L0 2 -co ro 2 ) 2 

+2e„ Ko 2 - ft) 2 XftV, 2 - ft> 2 ) 4 % (?> ft> +' r ) 

+(ax l0 2 -(0 2 )\{Anx K (q,(0+ir)} 2 

+{4%(«,fi)+iT)} 2 ]. (20) 

Section 3 below contains the procedures for extracting 
electron densities from the Fermi energies. 



3. Numerical Results - Electron Density 
and Fermi Energy 

Because the Fermi energy is one of the variables for 
calculating the Raman line shape, we give the numeri- 
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cal results for calculating the Fermi energy in terms of 
the BGN and PDOS models and the analytic expres- 
sions that relate the carrier concentrations to the Fermi 
energy. 

3.1 BGN Model 

We solve self-consistently by means of an iterative 
procedure, Eq. (3) with the distorted-perturbed carrier 
densities of states p Q {E) and py(E) from Ref. [9] used 
in the numerical integrations for n and h given by Eq. 
(1). The p Q {E) and py(E) are in tabular form and have 
both localized and continuum states. Their associated 
band edges contain the many body effects related to 
exchange and correlation. The independent variables 
are the temperature T and donor density 7V D . The Fermi 
energy is varied for a given temperature until Eq. (3) is 
satisfied to within an error of plus or minus 1 0~ 3 x N D . 

Numerical fitting procedures then give a closed-form 
analytic expression for the dependence of carrier con- 
centration /? BGN on the Fermi energy E F ; namely, 

1o &oKgn cm 3 ) 

= a BGN0 +a BGNl^F + a BGN2 ^F + — + °BGHl ^V '" ' \^) 

During the fitting analyses, we rely substantially on 
graphics and keep the number of fitting parameters to a 
minimum, subject to the constraint that the residual 
standard deviation 5 res is acceptably small; i.e., 5 res < 
0.01. The standard deviation is a measure of the "aver- 
age" error by which a fitted model represents a set of 
data points and thereby is a metric for assessing the 
quality of the fit. A smaller S res indicates a better fit. The 
residual standard deviation for a model Y f =f(Z) is 






— f 



= J y L(Y J -Yj) 2 /(N-P)], 



(22) 



where Y f are the calculated data values, the Yj are the 
predicted values from the fitted model, TV is the total 
number of data points (here N = 28), and P is the total 
number of parameters to be fitted in the model. We use 
the NIST-developed DATAPLOT [15] software for 
both the exploratory graphics and for the statistical 
analyses. 

In addition, we compare the above BGN model 
results with an equivalent two-band PDOS2 model for 
which there is no bandgap narrowing and for which the 
tabular carrier densities of states p c {E) and p v {E) from 
Ref. [9] are replaced with 



Pc(E)- 



N c 4kV^E 



(87T 3 )(r I '2m c m Q y 



and 



N 4jzV^[e 

Pv(^)= 3x^2/0 ^7T' res P ectivel Y- ( 23 ) 

(8tt )(H /2m Y m Q ) 

The corresponding polynomial fit for the carrier con- 
centration n versus Fermi energy E F relation is then 
denoted by 

log 10 (« cm 3 ) = a 00 + a 01 ^ + a^P^ +... + a Ql E l F ... . (24) 

Figure 1 compares the calculated carrier concentra- 
tion n BGN and n data for 28 values of donor densities 
between 10 16 cm -3 and 10 19 cm -3 . For a given Fermi 
energy, the electron densities predicted by the BGN 
model (solid-black curve) are typically factors of two 
smaller than the electron densities predicted by an 
equivalent parabolic two-band model (dashed-blue 
curve). 

3.2 PDOS Models 

We solve self-consistently, by means of an iterative 
procedure, Eq. (3) with Eq. (5). The independent vari- 
ables are the temperature T and donor density 7V D . 

We give here the results for fitting the logarithm to 
the base 10 of the total electron density and the electron 
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Figure 1. The calculated electron density «bgn f rom the BGN model 
as a function of Fermi energy for n-type GaAs at 300 K is given by 
the solid-black curve. The calculated electron density n from the 
two-band PDOS2 model, no bandgap narrowing model, as a function 
of Fermi energy for n-type GaAs at 300 K is given by the dashed- 
blue curve. The Fermi energy is relative to the majority conduction 
band edge at the r symmetry point in the first Brillouin zone. 



216 



Volume 112, Number 4, July-August 2007 

Journal of Research of the National Institute of Standards and Technology 



densities in each of the three conduction sub-bands at 
r, L, and X, n T , n L> and n X9 respectively, to polynomi- 
als in E F , namely 

log 10 O t cm 3 ) = a t0 +a tl E F +a a E F 2 +... + a tl E l F ..., (25) 

log 10 (/? r cm 3 ) = a ro + a n E F + a r2 E^ +...+a rl E l F ..., (26) 

log 10 (n L cm 3 )=a Ui +a u E F +a L2 E F +,„+a Ll E F ,„, and 

(27) 

log 10 (n x cm 3 ) =a xo +a xl E F + a X2 E F +.„+a xl E l F ... . (28) 

Figures 2 and 3 give the calculated electron densities 
as functions of the Fermi energy The corresponding fit- 
ted curves are not shown in Figs. 2 and 3 because the 
pairs of curves (calculated and fitted) for each of the 
electron densities n u n T , n h9 and n X9 lie on top of one 
another to within the line widths of each curve. Also, 
since the screening radii for the carriers from Eq. (2) 
are not needed when interpreting the proposed meas- 
urements considered here, the corresponding screening 
radii are not presented in this paper. Figure 2 compares 
the results from the PDOS2 model (dashed-blue curve 
with long spaces), PDOS2NPG model (solid-green 
curve), and PDOS4 model (dashed-red curve with short 
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Figure 2. Comparisons among three PDOS models with and with- 
out the non-parabolic factor £ in Eq. (5) for the electron energy dis- 
persion E cT (k). The solid-green curve shows the results for an equiv- 
alent two-band model that has one equivalent r conduction band 
with non-parabolicity factor £ and one equivalent r valence band; 
that is, the PDOS2NPG model. The dashed-blue curve with long 
spaces shows the results for an equivalent parabolic two-band model 
PDOS2 (£ = 0). The dashed-red curve short spaces shows the results 
for the four-band PDOS4 model. The Fermi energy is relative to the 
majority conduction band edge at the r symmetry point in the first 
Brillouin zone. 
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Figure 3. The calculated electron densities n t ,n r ,n L , and n x from 
the four-band PDOS4 model as functions of the Fermi energy for n- 
type GaAs at 300 K. The Fermi energy is relative to the majority con- 
duction band edge at the r symmetry point in the first Brillouin zone. 



spaces). The data for the blue curve in Fig. 1 is the same 
data for the blue curve in Fig. 2. Figure 3 shows the 
results for the four-band PDOS4 model. Unlike GaSb, 
most of the electrons are in the conduction sub-band at 
r . The conduction sub-band at L is only weakly popu- 
lated at the highest Fermi energies, and the conduction 
sub-band atXis negligible. 

3.3 Polynomial Fits - Closed-Form Analytic 
Expressions 

Tables 6 to 1 1 give the fitting parameters for polyno- 
mial fits to log 10 O BGN cm 3 ), log 10 (n cm 3 ), log 10 (« t cm 3 ), 
logioOr cm 3 ), log 10 (w L cm 3 ), and log 10 (n x cm 3 ) as shown, 
respectively, in Eqs. (21) and (24) to (28) and the cor- 
responding residual standard deviations *S res . In general, 
the values of S ms decrease monotonically with increas- 
ing number / of terms in these polynomials. But, care 
must be taken to avoid fitting noise in data sets. The 
general guideline for many data sets is that when the 
/-ratio (i.e., the absolute value of the ratio of the esti- 
mated parameter value divided by its estimated stan- 
dard deviation) is less than about 2, then the rate of 
decrease in S KS with increasing / tends to decrease, and 
proceeding with higher / values probably is not war- 
ranted. Because the changes in values of S rc$ between 
/ = 4 and 1=5 are not experimentally significant, we 
use the fitting parameters for the quartic case / = 4 in 
Eqs. (21) and (24). Similary, because the changes in 
values of S rcs between / = 3 and / = 4 are not experimen- 
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Table 6. Bandgap narrowing BGN model for log 10 (« BGN cm ). The five fitting parameters for a quar- 
tic polynomial fit Eq. (21) of the theoretical calculation for the equivalent conduction band electron 
density in n-type, zinc blende GaAs at 300 K as a function of the Fermi energy relative to the bottom 
of the equivalent conduction band at r. This quartic polynomial fit, which represents the theoretical 
results for Eq. (3), is valid only when -0.067 eV <E F < 0.286 eV. The *-ratio is the absolute value of 
the estimated fitting parameter \a BGN[ \ divided by its estimated standard deviation. The residual stan- 
dard deviation is S^ = 0.0130. 



Fitting parameter 


Estimated value 


Estimated standard 
deviation 


Units 


^-ratio 


a BGN0 


17.3292 


0.4191 x 10" 2 




4.135 x 10 3 


^BGNl 


13.1545 


0.7646 x 10" 1 


eV 1 


1.72 x 10 2 


a BGN2 


-37.4789 


1.473 


eV" 2 


25.45 


a BGN3 


26.5678 


12.86 


eV" 3 


2.067 


a BGN4 


53.7760 


29.43 


eV^ 


1.827 



Table 7. Two-band, no bandgap narrowing PDOS2 model for log 10 (« cm ). The five fitting parame- 
ters for a quartic polynomial fit Eq. (24) of the theoretical calculation for the L sub-band electron den- 
sity in n-type, zinc blende GaAs at 300 K as a function of the Fermi energy relative to the bottom of 
the conduction r sub-band. This quartic polynomial fit, which represents the theoretical results for Eq. 
(3), is valid only when -0.0974 eV <E F < 0.250 eV. The ^-ratio is the absolute value of the estimated 
fitting parameter \a Qi \ divided by the its estimated standard deviation. The residual standard deviation 
is S =0.0122. 



Fitting parameter 


Estimated value 


Estimated standard 
deviation 


Units 


^-ratio 


a oo 


17.5156 


0.3411 x 10" 2 




5.136 xlO 3 


a 0l 


12.8520 


0.8063 x 10 _1 


eV 1 


3.1 x 10 2 


a 02 


-34.5783 


0.6672 


eV" 2 


51.83 


a 03 


-27.3728 


9.166 


eV" 3 


2.986 


a m 


223.720 


29.54 


eV^ 


7.573 



Table 8. Four-band PDOS4 model for total electron density log 10 (n t cm ). The four fitting parameters 
for a cubic polynomial fit Eq. (25) of the theoretical calculation for the total electron density in n-type, 
zinc blende GaAs at 300 K as a function of the Fermi energy relative to the bottom of the conduction 
T sub-band. This cubic polynomial fit, which represents the theoretical results for Eq. (3), is valid only 
when -0.0974 eV < E F < 0.229 eV. The *-ratio is the absolute value of the estimated fitting parameter 
\a^\ divided by its estimated standard deviation. The residual standard deviation is S KS = 0.0234. 

Fitting parameter Estimated value Estimated standard Units ^-ratio 

deviation 

a, 



ti 

a l3 



17.5198 


0.6519 x 10" 2 




2.687 x 10 3 


12.3461 


0.9401 x 10" 1 


eV" 1 


1.31 x 10 2 


-35.5542 


1.322 


eV 2 


26.90 


40.1630 


6.488 


eV 3 


6.190 
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Table 9. Four-band PDOS4 model for electron density in the r sub-band log 10 (« r cm ). The four fit- 
ting parameters for a cubic polynomial fit Eq. (26) of the theoretical calculation for the r sub-band 
electron density in n-type, zinc blende GaAs at 300 K as a function of the Fermi energy relative to the 
bottom of the conduction r sub-band. This cubic polynomial fit, which represents the theoretical results 
for Eq. (3), is valid only when -0.0974 eV < E F < 0.229 eV. The /-ratio is the absolute value of the esti- 
mated fitting parameter \a Vl \ divided by its estimated standard deviation. The residual standard devia- 
tion is &= 0.0214. 



Fitting parameter 


Estimated value 


Estimated standard 
deviation 


Units 


/-ratio 


a ro 


17.5187 


0.5978 x 10" 2 




2.930 x 10 3 


a n 


12.3944 


0.8621 x 10" 1 


eV 1 


1.438 x 10 2 


a n 


-35.1795 


1.212 


eV" 2 


29.03 


«n 


33.7647 


5.950 


eV 3 


5.675 



14.2422 


0.3269 x 10" 3 




4.357 x 10 4 


16.8180 


0.4714 x 10" 2 


eV 1 


3.568 x 10 3 


0.162418 


0.6627 x 10" 1 


eV 2 


2.451 


-2.1988 


0.3253 


eV 3 


6.759 



Table 10. Four-band PDOS4 model for electron density in the L sub-band log 10 (w L cm )• Tne f° ur f" lt_ 
ting parameters for a cubic polynomial fit Eq. (27) of the theoretical calculation for the L sub-band 
electron density in n-type, zinc blende GaAs at 300 K as a function of the Fermi energy relative to the 
bottom of the conduction r sub-band. This cubic polynomial fit, which represents the theoretical results 
for Eq. (3), is valid only when -0.0974 eV < E F < 0.229 eV. The /-ratio is the absolute value of the esti- 
mated fitting parameter \a u \ divided by its estimated standard deviation. The residual standard devia- 
tion is S KS = 0.001 17. 

Fitting parameter Estimated value Estimated standard Units /-ratio 

deviation 

«L1 

a u. 

a l3 



Table 11. Four-band PDOS4 model for the electron density X sub-band log 10 (« x cm 3 ). The four fit- 
ting parameters for a cubic polynomial fit Eq. (28) of the theoretical calculation for the X sub-band 
electron density in n-type, zinc blende GaAs at 300 K as a function of the Fermi energy relative to the 
bottom of the conduction r sub-band. This cubic polynomial fit, which represents the theoretical results 
for Eq. (3), is valid only when -0.0974 eV < E F < 0.229 eV. The /-ratio is the absolute value of the esti- 
mated fitting parameter {a^ divided by its estimated standard deviation. The residual standard devia- 
tion is S res = 0.000 019 5. 

Fitting parameter Estimated value Estimated standard Units /-ratio 

deviation 

a xo 
a xl 

a X2 
a X3 



11.2919 


0.5452 x 10" 5 




2.071 x 10 6 


16.7999 


0.7862 x 10"> 


eV" 1 


2.137 x 10 5 


-0.542 16 x 10" 3 


0.1105 xlO -2 


eV 2 


0.4905 


-0.4041 x 10" 2 


0.5426 x 10" 2 


eV 3 


0.7448 
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tally significant, we use the fitting parameters for the 
cubic 1 = 3 case in Eqs. (25) to (28). 



4. Conclusions 

The foregoing theory for extracting electron densi- 
ties from Raman spectra is unique in two ways: 1) It 
treats the many-body effects self-consistently and 2) it 
is valid at room temperature for arbitrary values of the 
ratio R = (Q 2 /oc). When high concentrations of carriers 
exists, this theory and its associated numerical proce- 
dures for determining carrier concentrations from 
Fermi energies are necessary for interpreting room tem- 
perature Raman spectra self-consistently for arbitrary 
values of the ratio R, even for R ~ 1. The BGN models 
presented here include many-body quantum effects and 
bandgap narrowing. However, obtaining Fermi ener- 
gies from experimental Raman spectra is beyond the 
scope of the present paper and involves the following 
steps to implement contactless measurements of carrier 
densities: 

1) Numerical evaluation of Eq. (17) with adaptive 
grids to treat the integrable singularities of the inte- 
grands, 

2) Three-dimensional visualization of computer results 
for predicted Raman line shapes as functions of Fermi 
energy and frequency, and 

3) Iteration of the predicted Raman line shape with the 
Fermi energy as a variation parameter to give the best 
self-consistent fit to the measured Raman line shape. 
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